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Abstract 

Results are presented from a study of the structure of high energy hadronic 
events recorded by the L3 detector at yfs > 192 GeV. The distributions of sev- 
eral event shape variables are compared to resummed 0(a^) QCD calculations. 
We determine the strong coupling constant at three average centre-of-mass ener- 
gies: 194.4, 200.2 and 206.2 GeV. These measurements, combined with previous 
L3 measurements at lower energies, demonstrate the running of a s as expected in 
QCD and yield a s (m z ) = 0.1227 ± 0.0012 ± 0.0058, where the first uncertainty is 
experimental and the second is theoretical. 
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1 Introduction 



The measurement of the energy dependence of the strong coupling constant, a s , is an important 
test of Quantum Chromo-Dynamics (QCD). Hadronic events produced in e + e~ annihilation 
offer a clean environment to perform such measurements. The high energy phase of the LEP 
collider gives a unique opportunity to measure QCD observables over a wide energy range 
and to perform a precise test of the energy dependence of the strong coupling constant. In 
addition, the study of hadronic events allows to check the validity of the QCD models used for 
background modelling in other studies, such as new particle searches. 

In its last two years, the LEP collider operated at various centre-of-mass energies, y/s, 
between 192 and 208 GeV. These are grouped in three samples of average centre-of-mass 
energies 194.4, 200.2 and 206.2 GeV, corresponding to the ranges 192 < y/s < 196 GeV, 
200 < y/7s < 202 GeV and ^/s~ > 202 GeV respectively. 

We report on measurements of five event shape distributions using 436.8 pb _1 of data 
collected with the L3 detector [1] at these centre-of-mass energies, as detailed in Table 1. To 
allow a direct comparison with our earlier QCD studies at lower energies [2-6], we follow an 
identical analysis procedure. The value of a s is extracted in each energy range by comparing 
the measured event shape distributions with predictions of second order QCD calculations [7] 
supplemented by resummed leading and next-to- leading order terms [8-11]. These values are 
used, together with previous L3 measurements at lower effective centre-of-mass energies, from 
30 GeV to 189 GeV, to study the energy evolution of a s . 

2 Event Selection 

The criteria for the selection of e + e~ — > qq — > hadrons events are identical to those used in 
our previous QCD study at y/s =189 GeV [6]. They are based on the measured total visible 
energy, E vis , the energy imbalances parallel, E\\, and perpendicular, E±, to the beam direction 
and on the cluster multiplicity, N d . These variables are constructed using energy clusters in 
the electromagnetic and hadronic calorimeters with a minimum energy of 100 MeV. 

The efficiency of the selection criteria and the purity of the data sample are estimated 
using Monte Carlo events for the process e + e~ — > qq(7) generated by the KK2F [12] program, 
interfaced with JETSET PS [13] routines to describe the QCD parton shower evolution and 
hadronisation. The events are then passed through the L3 detector simulation [14]. The 
KK2F generator is chosen for its improved simulation of the initial state radiation (ISR) as 
compared to the PYTHIA [13] model previously used. Background events are simulated with 
PYTHIA for two-photon events and Z-pair production, KORALZ [15] for the r + T'( / ~f) final 
state, BHAGENE [16] and BHWIDE [17] for Bhabha events and KORALW [18] for W-pair 
production. 

Hadronic events with hard ISR photons, where the mass of the hadronic system is close 
to roz, are considered as background if the photon energy exceeds 0.18\/s. This important 
background is reduced to less than 8% of the selected events by applying a cut in the plane 
of \E\\\/E vis vs. E vis /y/s. Additional background arises from W- and Z-pair production. A 
substantial fraction of these events is removed by a specific selection [5] that forces a 4- jet 
topology using the Durham algorithm [19], and applies cuts on the cluster multiplicity, N d > 
40, the jet resolution parameter, y^ 4 >0.0025, the energy of the least energetic jet and on the 
energy fraction carried by the two most energetic jets. The cuts are optimised to maximise the 
product of efficiency and purity at each energy point. After selection, the W-pair background 
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amounts to 6.4% at y/s = 192 GeV and increases to 10.3% at y/s = 208 GeV. The Z-pair 
background is below 0.8%. The selection efficiency, purity and number of selected events for 
the three energy points are summarised in Table 1. 

3 Measurement of Event Shape Variables 

The measured global event shape variables are the thrust [20], T, the scaled heavy jet mass [21], 
p, the total, B T) and wide, B W) jet broadening variables [9] and the C-parameter [22]. The 
first four observables are defined in terms of the particle four-momenta, while the C-parameter 
is derived from the spherocity tensor: 
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where the sums run over all particles and p a is the momentum vector of the particle a. The 
C-parameter is defined in terms of the eigenvalues, Ai, A2 and A3, of y , as: 

C = 3(A 1 A 2 + A 2 A3 + A 3 A 1 ) . 

For all five variables, improved analytical QCD calculations [8-11] are available. The calcu- 
lations used here for the jet broadening variables [11] are improved as compared to the previous 
predictions [9] by a better treatment of quark recoil effects. 

After background subtraction, the measured distributions are corrected bin-by-bin for de- 
tector effects, acceptance and resolution. The correction factors are the ratios of Monte Carlo 
distributions at detector level to the distributions at particle level which include all stable 
charged and neutral particles. 1 - ) The data are also corrected bin-by-bin for initial and fi- 
nal state photon radiation using Monte Carlo distributions at particle level with and without 
radiation. 

Figures 1 and 2 show the thrust and wide jet broadening distributions corrected to the 
particle level. The data are compared with the JETSET, HERWIG [23] and ARIADNE [24] 
QCD models. These models, based on an improved leading-logarithmic approximation parton 
shower, including QCD coherence effects, are tuned to reproduce the global event shape dis- 
tributions and the charged particle multiplicity distribution measured at 91.2 GeV [25]. At 
and above centre-of-mass energies of 200.2 GeV some discrepancies appear for specific values 
of the observables. Several studies are performed to investigate the reason for such effects. 
The observed structures in the global event shape distributions are found to depend neither on 
time nor on detector geomety. The effects of these discrepancies are taken into account in the 
determination of the systematic uncertainties. 

The two main sources of systematic uncertainty in the event shape variable distributions 
are those on detector correction and background estimation. These uncertainties are estimated 
by repeating the measurement with different analysis criteria and correction procedures [5,6]. 
The uncertainty in the detector correction is estimated with the following tests: 

• The effect of different particle fluxes in correcting the measured distribution is estimated 
by using the HERWIG Monte Carlo program instead of JETSET to simulate the signal. 
Half of the difference obtained with these two models is taken as a systematic uncertainty. 

^All weakly decaying light particles with mean lifetime larger than 3.3 x 10~ 10 s are considered stable. 
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• The definition of reconstructed objects used to calculate the observables is changed from 
calorimetric clusters to a non-linear combination of energies of charged tracks and calori- 
metric clusters. 

• The acceptance is reduced by restricting the analysis to events in the central part of the 
detector, |cos(0t)| < 0.7, where 9t is the polar angle of the thrust axis relative to the 
beam direction. In this region a better energy resolution is found. 

Half of the maximum spread between the latter two tests and the original analysis is assigned 
as a systematic uncertainty. 

The uncertainty on the background composition of the selected sample is estimated by 
repeating the analysis with: 

• an alternative criterion to reject the hard initial state photon events based on a cut on 
the kinematically reconstructed effective centre-of-mass energy {\fs i < 0.92); 

• a variation of the W + W~ background estimate by applying a full subtraction of the 
W-pair contribution without preliminary event rejection; 

• a variation of the estimated two-photon interaction background by ± 30%. The Monte 
Carlo program used to model two-photon interactions is also changed from PYTHIA to 
PHOJET [26]. 

For the first two studies, half of the difference between the results of the original analysis and 
of the systematic check is taken as the systematic uncertainty. In the two-photon case, half 
of the maximum spread between the new results and the original analysis is considered as a 
systematic uncertainty. The statistical component of each systematic uncertainty is estimated 
and removed following the procedure described in reference [6]. The systematic uncertainties 
obtained from the different sources are then combined in quadrature. For y/s < 196 GeV, 
the uncertainties due to the backgrounds are the most important ones. They are 2 — 3 times 
larger than the uncertainties due to detector corrections. For y/s > 200 GeV, the uncertainty 
in the detector correction gives the largest systematic contribution, dominated by the effect 
of reducing the event thrust acceptance in the central part of the detector, but decreases for 
< v/s >= 206.2 GeV. 

An important test of QCD models is the comparison of the energy evolution of the means 
of the event shape variables. The mean values of the five variables obtained at y/s > 192 GeV 
are given in Table 2. Figure 3 shows the evolution of < 1 — T > and < By/ > as a function 
of y/s. Also shown are the energy dependences of these quantities as predicted by JETSET 
PS, HERWIG, ARIADNE, COJETS [27] and JETSET ME, with an 0{a 2 s ) matrix element 
implementation. For high energies, the JETSET ME and COJETS models are not favoured by 
the data. 

4 Determination of ct s 

The QCD predictions for the five event shape observables are based on O(a^) perturbative 
QCD calculations with resummed leading and next-to-leading order terms. To compare these 
calculations at parton level with the experimental distributions, the effects of hadronisation 
and decays are corrected for with a folding matrix [4] calculated using the JETSET PS Monte 
Carlo program. 
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To determine a s at each energy point, the measured distributions are fitted in the ranges 
given in Table 3 to the analytical predictions, using the modified-log(R) matching scheme [10] 
after corrections for hadronisation effects. Figure 4 shows the experimental data together with 
the result of the QCD fits for the five variables at < \/s > = 206.2 GeV. 

The a s measurements at the three energy points are summarised in Table 3 together with 
their experimental and theoretical uncertainties. The former includes the statistical and the 
experimental systematic uncertainties discussed above. The latter is obtained from estimates [4] 
of the hadronisation uncertainty and of the uncalculated higher orders in the QCD predictions. 

The hadronisation uncertainty is obtained from the variation in the fitted value of a s due to 
hadronisation corrections determined by comparing JETSET with HERWIG and ARIADNE 
models and changing the JETSET fragmentation parameters, b, a q and A LLA within their 
errors [25] and turning off Bose-Einstein correlations. The most important variation comes from 
the change in the fragmentation models and is taken as an estimate of the overall hadronisation 
uncertainty. 

The uncertainty coming from uncalculated higher orders in the QCD predictions is esti- 
mated in two independent ways: by varying the renormalisation scale, //, and by changing the 
matching scheme. The scale uncertainty is obtained by repeating the fit for different values of 
the renormalisation scale in the interval 0.5^ < // < 2y/s. The matching scheme uncertainty 
is obtained from half of the maximum spread given by different algorithms [10]. The largest of 
these uncertainties is assigned as the theoretical uncertainty due to uncalculated higher orders. 

To obtain a combined value for the strong coupling constant, we take the unweighted av- 
erage of the five a s values. The overall theoretical uncertainty is obtained from the average 
hadronisation uncertainty added in quadrature to the average higher order uncertainty. A 
cross-check of this theoretical uncertainty is obtained from a comparison of a s measurements 
from the various event shape variables which are expected to be differently affected by higher 
order corrections and hadronisation effects. Half of the maximum spread in the five a s values 
is found to be consistent with the estimated theoretical uncertainty. 

Earlier L3 measurements at ^/s = mz and at reduced centre-of-mass energies determined 
a s from four event shape variables only: T, p, Bt and Bw, the resummed calculation for the 
C-parameter not being available. We have determined a s at these lower energies from the 
C-parameter and the values are now included in the overall mean a s and listed in Table 4. 

The improved theoretical predictions for the jet broadening variables are used to update 
our previously published a s results at effective centre-of-mass energies from 30 GeV up to 189 
GeV, as listed in Table 4. The mean a s values from the five event shape distributions are given 
in Table 5 together with the experimental and theoretical uncertainties. 

Figure 5a compares the energy dependence of the measured a s values with the prediction 
from QCD. The theoretical uncertainties are strongly correlated between these measurements. 
Hence, the energy dependence of a s is investigated using only experimental uncertainties. The 
experimental systematic uncertainties on a s are partially correlated. The background uncer- 
tainties are similar for data points in the same energy range but differ between the low energy, Z 
peak and high energy data sets. The sixteen measurements in Figure 5a are shown with exper- 
imental uncertainties only, together with a fit to the QCD evolution equation [28] with a s (mz) 
as a free parameter, that takes into account the correlation between the various measurements. 
The covariance matrix for the fit is obtained as follows: 

• The statistical uncertainties are taken as uncorrelated. 

• The experimental systematic uncertainties are assumed to be uncorrelated between the 
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three data sets and to have a minimum overlap correlation between different energies 
within the same data set. This definition consists of assigning to the covariance matrix 
element the smallest of the two squared uncertainties. 

The fit gives a x 2 °f 17.9 for 15 degrees of freedom corresponding to a confidence level of 
0.27 yielding a value of a s : 

a s (m z ) = 0.1227 ± 0.0012 ± 0.0058 . 

The first uncertainty is experimental and the second theoretical. The latter is obtained 
from the result of a fit which includes the theoretical uncertainties and their correlations. The 
covariance matrix is here defined assuming a minimum overlap correlation between energies of 
the hadronisation as well as the uncalculated higher order uncertainties. The hadronisation 
uncertainty contribution to the total theoretical uncertainty is ±0.0026. 

A fit with constant a s gives a \ 2 °f 51.7 for 15 degrees of freedom. These measurements 
support the energy evolution of the strong coupling constant predicted by QCD. The apparent 
increase of the a s values obtained at > 194 GeV compared to the QCD evolution curve is 
related to the structures seen in the event shape distributions discussed above. 

Figure 5b summarises the a s values determined by L3 from the measurement of the r branch- 
ing fractions into leptons [29], Z lineshape [30] and event shape distributions at various energies, 
together with the QCD prediction obtained from the fit to the event shape measurements only. 
The band width corresponds to the evolved uncertainty on « s (mz). All the measurements are 
consistent with the energy evolution of the strong coupling constant predicted by QCD. The 
uncertainties on these measurements are dominated by the theoretical uncertainty coming from 
the unknown higher order contributions in the calculations. An improved determination of a s 
from these measurements thus awaits improved theoretical calculations of these corrections. 
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< y/s > (GeV) 


194.4 


200.2 


206.2 


Integrated Luminosity (pb _1 ) 


112.2 


117.0 


207.6 


Selection Efficiency (%) 


82.8 


85.7 


86.0 


Sample Purity (%) 


81.4 


80.6 


78.8 


Selected Events 


2403 


2456 


4146 



Table 1: Summary of integrated luminosity, selection efficiency, sample purity and number of 
selected hadronic events at the average centre-of-mass energies used in this analysis. 



< Vs > 


194.4 GeV 


200.2 GeV 


206.2 GeV 


< (1 -T) > 


0.0551 ± 0.0021 ± 0.0009 


0.0582 ± 0.0021 ± 0.0015 


0.0569 ± 0.0017 ± 0.0016 


<P> 


0.0439 ± 0.0014 ± 0.0007 


0.0464 ± 0.0014 ± 0.0015 


0.0455 ± 0.0011 ± 0.0011 


< B T > 


0.0920 ± 0.0022 ± 0.0023 


0.0950 ± 0.0021 ± 0.0025 


0.0938 ± 0.0017 ± 0.0015 


< B w > 


0.0663 ± 0.0014 ± 0.0009 


0.0688 ± 0.0013 ± 0.0016 


0.0682 ± 0.0011 ± 0.0009 


<C> 


0.2158 ± 0.0058 ± 0.0023 


0.2244 ± 0.0059 ± 0.0068 


0.2195 ± 0.0049 ± 0.0035 



Table 2: Mean values of the five event shape variables at different energy points. The first 
uncertainty is statistical and the second systematic. 
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1-T 


P 


Bj' 


B w 


C 


Fit Range 


0.00-0.25 


0.00-0.20 


0.02-0.26 


0.015-0.210 


0.05-0.50 


a s (194.4 GeV) 


0.1168 


0.1096 


0.1152 


0.1071 


0.1130 


Statistical uncertainty 


±0.0019 


±0.0017 


±0.0015 


±0.0017 


±0.0023 


Systematic uncertainty 


±0.0015 


±0.0014 


±0.0015 


±0.0013 


±0.0024 


Overall experimental uncertainty 


±0.0024 


±0.0022 


±0.0021 


±0.0021 


±0.0033 


Overall theoretical uncertainty 


±0.0056 


±0.0039 


±0.0065 


±0.0062 


±0.0056 


xVd.o.f. 


2.2 / 9 


10.1 / 13 


20.2 / 11 


9.7 / 12 


3.8 / 8 


a s (200.2 GeV) 


0.1178 


0.1114 


0.1164 


0.1088 


0.1147 


Statistical uncertainty 


±0.0019 


±0.0017 


±0.0015 


±0.0017 


±0.0024 


Systematic uncertainty 


±0.0027 


±0.0028 


±0.0018 


±0.0014 


±0.0016 


Overall experimental uncertainty 


±0.0033 


±0.0033 


±0.0023 


±0.0022 


±0.0029 


Overall theoretical uncertainty 


±0.0059 


±0.0034 


±0.0062 


±0.0062 


±0.0057 


xVd.o.f. 


7.3 / 9 


6.8 / 11 


9.6 / 11 


10.3 / 12 


2.9 / 8 


a s (206.2 GeV) 


0.1173 


0.1119 


0.1163 


0.1077 


0.1130 


Statistical uncertainty 


±0.0014 


±0.0013 


±0.0012 


±0.0014 


±0.0019 


Systematic uncertainty 


±0.0016 


±0.0014 


±0.0017 


±0.0013 


±0.0021 


Overall experimental uncertainty 


±0.0021 


±0.0019 


±0.0021 


±0.0019 


±0.0028 


Overall theoretical uncertainty 


±0.0057 


±0.0034 


±0.0065 


±0.0062 


±0.0053 


xVd.o.f. 


7.5 / 9 


7.7 / 13 


5.9 / 11 


7.8 / 12 


3.4 / 8 



Table 3: Values of a s measured at < \/s > = 194.4, 200.2 and 206.2 GeV from fits of the event 
shape variables. The fit ranges, the estimated experimental and theoretical uncertainties and 
the x 2 /d.o.f. of the fit are also given. 



< Vs > 
(GeV) 


a s measurement 


from Bt 


from B w 


from C 


41.4 





1401±0 


0063±0 


0119 





1380±0.0067±0 


0091 





1371±0 


0070±0 


0102 


55.3 





1321±0 


0070±0 


0099 





1191±0.0072±0 


0088 





1197±0 


0086±0 


0118 


65.4 





1354±0 


0067±0 


0106 





1190±0.0062±0 


0086 





1258±0 


0039±0 


0108 


75.7 





1296±0 


0074±0 


0097 





1068±0.0065±0 


0084 





1143±0 


0072±0 


0094 


82.3 





1270±0 


0079±0 


0095 





1083±0.0067±0 


0087 





1153±0 


0060±0 


0091 


85.1 





1259±0 


0069±0 


0095 





1092±0.0080±0 


0091 





1115±0 


0045±0 


0089 


91.2 





1222±0 


0020±0 


0080 





1196±0.0022±0 


0052 





1170±0 


0016±0 


0076 


130.1 





1178±0 


0033±0 


0064 





1089±0.0031±0 


0088 





1151±0 


0040±0 


0066 


136.1 





1166±0 


0035±0 


0064 





1072±0.0041±0 


0078 





1089±0 


0047±0 


0076 


161.3 





1123±0 


0042±0 


0067 





1058±0.0059±0 


0068 





1043±0 


0060±0 


0057 


172.3 





1105±0 


0063±0 


0061 





1062±0.0050±0 


0065 





1121±0 


0068±0 


0057 


182.8 





1145±0 


0022±0 


0060 





1045±0.0021±0 


0071 





1081±0 


0029±0 


0054 


188.6 





1153±0 


0018±0 


0067 





1063±0.0017±0 


0078 





1118±0 


0023±0 


0055 



Table 4: Updated a s measurements from the jet broadening distributions and the C-parameter 
measured at y/s < 189 GeV. The first uncertainty is experimental and the second theoretical. 
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< ^fs > 




a s 


measurement from T, p, 


B T , B 


w, C 


(GeV) 


a s 


stat 


syst 


hadr. 


hi. 


order 


41.4 





1418 


±0.0053 


±0.0030 


±0 


0055 


±0 


0085 


55.3 





1260 


±0.0047 


±0.0056 


±0 


0066 


±0 


0062 


65.4 





1331 


±0.0032 


±0.0042 


±0 


0059 


±0 


0064 


75.7 





1204 


±0.0024 


±0.0059 


±0 


0060 


±0 


0053 


82.3 





1184 


±0.0028 


±0.0053 


±0 


0060 


±0 


0051 


85.1 





1152 


±0.0037 


±0.0051 


±0 


0060 


±0 


0055 


91.2 





1210 


±0.0008 


±0.0017 


±0 


0040 


±0 


0052 


130.1 





1138 


±0.0033 


±0.0021 


±0 


0031 


±0 


0046 


136.1 





1121 


±0.0039 


±0.0019 


±0 


0038 


±0 


0045 


161.3 





1051 


±0.0048 


±0.0026 


±0 


0026 


±0 


0044 


172.3 





1099 


±0.0052 


±0.0026 


±0 


0024 


±0 


0048 


182.8 





1096 


±0.0022 


±0.0010 


±0 


0023 


±0 


0044 


188.6 





1122 


±0.0014 


±0.0012 


±0 


0022 


±0 


0045 


194.4 





1123 


±0.0018 


±0.0016 


±0 


0020 


±0 


0047 


200.2 





1138 


±0.0018 


±0.0021 


±0 


0020 


±0 


0046 


206.2 





1132 


±0.0014 


±0.0016 


±0 


0019 


±0 


0047 



Table 5: Combined a s values from the five event shape variables with their uncertainties. 
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Figure 1: Corrected distributions for the thrust, T, at < y/s >= 194.4, 200.2 and 206.2 GeV 
compared with QCD model predictions. The uncertainties shown are statistical only. 
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Figure 2: Corrected distributions for the wide jet broadening, Bw, at < y/s >= 194.4, 200.2 
and 206.2 GeV compared with QCD model predictions. The uncertainties shown are statistical 
only. 
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Figure 3: Distribution of a) < 1 — T > and b) < B w > as a function of the centre-of-mass 
energy, compared to several QCD models. Lower energy data [2-6] are also presented. The 
error bars include experimental systematic uncertainties. 
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Figure 4: Measured distributions of thrust, T, scaled heavy jet mass, p, total, B?, and wide, 
B w , jet broadenings, and C-parameter in comparison with QCD predictions at < y/s >=206.2 
GeV. The error bars include experimental systematic uncertainties. 
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Figure 5: Values of a s determined as a function of yt~s: a) from event shape distributions with 
experimental uncertainties only. The solid and dashed lines are fits with the energy dependence 
of a s as expected from QCD and with constant a s , respectively; b) from the measurement of 
the r branching fractions into leptons [29] , Z lineshape [30] and event shape distributions. The 
dashed line is a fit to the QCD evolution function to the measurements made from event shape 
variables. The band width corresponds to the evolved uncertainty on a s (m z ). 
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